Astronomy Sc Astrophysics manuscript no. 14612 © ESO 2010 

May 14, 2010 



d . 6 



Letter to the Editor 

Hier ist wahrhaftig ein Loch im Himmel ^ ^^ 

The NGC 1999 dark globule is not a globule 

T. Stanke\ A. M. Stutz^'^ J. J. Tobin^, B. Ali^ S. T. Megeath^, O. Krause^, H. Linz^, L. Allen^ E. Bergin^, N. Calvet"^, 

J. Di Francesco^'^, W. J. Fischer^, E. Furlan^^, L. Hartmann"^, T. Henning^, P. Manoj^^, S. Maret^^, J. Muzerolle^^, 

P. C. Myers^^, D. Neufeld^^ M. Osorio^^, K. Pontoppidan^^, C. A. Poteet^, D. M. Watson^^ and T. Wilson^ 

^ ESO, Karl-Schwarzschild-Strasse 2, 85748 Garching bei München, Germany; e-mail: tstanke@eso . org 
^^ ' ^ Max-Planck-Institut für Astronomie, Königstuhl 17, D-69117 Heidelberg, Germany 

^v^ \ ^ Department of Astronomy and Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA 

^ Department of Astronomy, University of Michigan, Ann Arbor, MI 48109, USA 
^. ^ NASA Herschel Science Center, CaUfornia Institute of Technology, 770 South Wilson Ave, Pasadena, CA91 125, USA 

Department of Physics and Astronomy, University of Toledo, 2801 West Bancroft Street, Toledo, OH 43606, USA 
National Optical Astronomy Observatory, 950 N. Cherry Ave., Tucson, AZ 85719, USA 

Department of Physics and Astronomy, University of Victoria, PO. Box 355, STN CSC, Victoria BC, V8W 3P6, Canada 
CN| . ^ National Research Council Canada, Herzberg Institute of Astrophysics, 5071 West Saanich Road, Victoria BC, V9E 2E7, Canada 

" JPL, California Institute of Technology, Mail Stop 264767, 4800 Oak Grove Drive, Pasadena, CA 91 109, USA 
Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627, USA 
^ , Laboratoire d'Astrophysique de Grenoble, Universite Joseph Fourier, CNRS, UMR 571, BP 53, F-38041 Grenoble, France 

^ . ^^ Space Telescope Science Institute, 3700 San Martin Dr., Baltimore, MD 21218, USA 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 

Department of Physics and Astronomy, Johns Hopkins University, 3400 North Charles Street, Baltimore, MD 21218, USA 
r"| , ^6 Instituto de Astrofisica de Andalucia, CSIC, Camino Bajo de Huetor 50, E- 18008, Granada, Spain 

rS |. ^'^ Division of Geological and Planetary Sciences 150-21, California Institute of Technology, Pasadena, CA 91 125, USA 

ö 

^ ■ ABSTRACT 

c/) . 

Cu . The NGC 1999 reflection nebula features a dark patch with a size of ~ 10,000 AU, which has been interpreted as a small, dense 

foreground globule and possible site of imminent star formation. We present Herschel PACS far-infrared 70 and 160yL/m maps, which 
reveal a flux deficit at the location of the globule. We estimate the globule mass needed to produce such an absorption feature to 
be a few tenths to a few M©. Inspired by this Herschel Observation, we obtained APEX LABOCA and SABOCA submillimeter 
continuum maps, and Magellan PANIC near-infrared images of the region. We do not detect a submillimer source at the location of 
the Herschel flux decrement; furthermore our observations place an upper limit on the mass of the globule of ~2.4-10~^ M©. Indeed, 
the submillimeter maps appear to show a flux depression as well. Furthermore, the near-infrared images detect faint background stars 
that are less aff'ected by extinction inside the dark patch than in its surroundings. We suggest that the dark patch is in fact a hole or 
cavity in the material producing the NGC 1999 reflection nebula, excavated by protostellar jets from the V 380 Ori multiple System. 
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1 . Introduction dust content absorbs the light of background stars. Dark clouds, 

ranging in mass and size from giant molecular clouds (GMCs) 

In the year 1774 C.E. Sir Friedrich Wilhelm Herschel first no- jo tiny globules, are the sites of star formation in our Galaxy. 

^ ticedpatchesofskyinthecoiistellationScorpiothatweredevoid NGC1999 is a small reflection nebula in the LDN1641 

Cd of Stars. Unable to find even the faintest star in these regions, his ^^^^^^ ^^ ^^^ Orion A GMC, located in a small group of 22 

sister Caroine reported htm to exclaim: Hier ist wahrhaftig ein _^^^^ ^equence stars and protostars, including the driving 

Loch im Himmel! ( Truly there is a hole in the sky here! ). ^^^^^^ ^f the prototypical Herbig-Haro objects HH 1 and HH 2 

These da rk areas are now kn own to be due to obscunng mate- (Megeath et al., in prep.). It is illuminated by the Herbig Ae/Be 

rial (e.g.,| Bamardetal. || l927|) : clouds of molecular gas, whose ^j^ v 380 Ori, a multiple System with a circumsystem disk 

where the primar y is a 100 Lg B9 st ar exhibiting strong emis- 

* Herschel is an ESA space observatory with science Instruments sion lines (e.g.. iLeinert et all [TÖÖtL HillenbrandeFüJ [Uli 

provided by European-led Principal Investigator consortia and with im- |Aiecian et alJ |2009). NGC 1999 features a compact (20-30") 

portant participation trom NASA . . ^ i / t- m i • i • j -i j • ^i t^ 

^.^ rr., 1 V- • • 1 1 1 -1 • , , A dark patch (see Fig. IIJ, which is descnbed in the literature as 

This publication mcludes data acquired with the Atacama j i i i i i . .• i •. r ^ r ^. / ITT — 7"n 

Pathfinder Experiment (APEX; proposal E-082.F-9807 and E-284.C- i ^il^l&l^, ^"^ "^T }. ". I . no^ ^u u u u ^^•^■' ^^ 

5015). APEX is a coUaboration between the Max-Planck-Institut für 1946, 1960|;| Warren-Smith et al. | 1980|) , but which has neverbeen 

Radioastronomie, the European Southern Observatory, and the Onsala stuaiea in greater detail berore. 

Space Observatory This paperincludes data gathered with the 6.5 meter We report here the serendipitous and surprising detection 

Magellan Telescopes located atLas Campanas Observatory, Chile. with PACS of a 70 and 160yum flux decrement at the location 
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Fig.l. NGC 1999 and HH 1/2 region (DSS, left) and HST 
F450W/F555W/F675W true color Image of the NGC 1999 dark 
patch. 




Fig. 2. NGC 1999/HH1/2 region as seen with Herschel PACS 
at 160yum (left, the box marks the region shown in Fig.O, and 
APEX SABOCA at 350yum (right, smoothed to 10'' resolution), 
all displayed with a logarithmic color Stretch. 



of NGC 1999. The maps show a dark patch against the nebu- 
lous far-IR emission that closely resembles the dark patch seen 
at visible wavelengths. The MSX, ISO, and Spitzer space tele- 
scopes have detected clouds in absorption at mid-IR (5-30/im) 
wave lengths against the diffuse IR background of the Galax y 
(e,g., 'Bacmann etäDl2000|; IStutz et al.1 [20081: iTobmet al.''2ÖT0; 
Stutz e t al. 2009^, in a few cases even at 70 //m f Stutz et al. 
2009bh . However, the detection of the NGC 1999 globule by 
PACS would be the first detection of a dark globule in absorp- 
tion at 160 /im. At wavelengths > 160yum, cold globules are nor- 
mally detected through the emission from cold (10-30 K) dust. 
Furthermore, absorption at these wavelengths requires extinc- 
tions in excess of lOOAy. 

Motivated by the possible discovery of a 160yum dark cloud, 
we obtained follow-up ground based submillimeter and near-IR 
(extinction) observations towards the globule. We did not de- 
tect the column density or mass of cold dust necessary to pro- 
duce such an absorption feature in the far-IR. We therefore con- 
clude that the dark globule is actually a cavity in the NGC 1999 
cloud carved by outflows from the V 380 Ori System. Thus, the 
Herschel telescope has discovered what truly is a hole in the sky. 

2. Observations and processing 

2.1. Herschel PACS observations 



NGC 1999 was observed with Herschel (|Pilbratt et al.1 l2010l) 
with the PACS Instrument (iPoglitsch et al.1 l2010l) at 70 
and 160yum on 2009 October 9 (OBSIDs 1342185551 and 
1342185552) as part of the Science Demonstration Phase ob- 
servations of the Herschel Orion Protostar Survey (HOPS; PI. 
S. T. Megeath). The data are presented in Figs.[2]and[3](see also 
iFischer et al.l l2010b . They were processed with the Herschel 
Common Software System (HCS© version 3.0 build 919, fol- 
lowing Standard Steps, using v. 4 of the calibration file. After ini- 
tial processing the two orthogonal scans were combined into a 
Single Observation. We used the same HCSS Software modules as 
the automatic PACS pipeline. We used two different approaches 
for combining the data into final mosaics: 
Methodl: Simple (''naive'') mapping with baseline removal. To 



preserve the extended emission we estimated the s ky and the cor- 
related signal drifts observed in PACS readouts (" Sauvage et al.l 
2010). To determine the baselines, we removed pixel-to-pixel 
electronic offsets using the median of the entire time stream, 
then calculated the median for each frame/readout. The median 
values were divided into 500 readout bins, and we took the min- 
imum value for each bin. This approach preserves the extended 
emission at all spatial scales, but does not remove 1/f noise. The 
final mosaic was created using the "photProject" routine. 
Methodl: Optimal mapping with MADmapß to remove the 
bolometer 1/f noise from the pixel time-line. The initial pro- 
cessing is the same as for Method 1. However, instead of 'phot- 
Project', MADmap is used to determine an optimal Solution for 
each sky pixel. This method preserves extended emission and 
removes the 1/f noise, at the price of a reduced sensitivity. 

2.2. Otherdata 

APEX SABOCA and LABOCA — We obtained submillimeter 
continuum maps using LABO CA and SABOCA on APEX. 
LABOCA (ISiringoetal.ll2009h is a -250 bolometer array op- 
erating at 870/im, with a spatial resolution of 19''. Observations 
were done on 2009 November 29 and 30 in fair conditions, 
with precipitable water vapour (PWV) values around 1mm 
(Tzenith,870 ~ 0.26). Wc uscd 3. combination of spiral and straight 
on-the-fly scans in order to recover extended emission. Data 
reduction was done with the BOA Software (Schuller et al., in 
prep.) following Standard procedures, including iterative source 
modeling. SABOCA (Siringo et al. 2010) is a 37 bolometer ar- 
ray operating at 350 yum, with a resolution of 7V8. Data were 
taken on 2009 December Ist and 3rd (NGC 1999) and on 2008 
October 7 and 8 (HH 1/2 area) with PWV<0.5 mm (Tzenith,350 < 
1). The observing and data reduction procedures were similar to 
those used for LABOCA. The 350 yum map is presented in Fig.O 
and both maps in Fig. lA.lK available in electronic form only). 
Magellan PANIC — NGC 1999 was observed with the PANIC 
near-IR camera (2' field of view) at Magellan on 2009 
December 4 in Photometrie conditions through //, K^ and H2 fil- 
ters in 0'.'3 to 0'.'4 seeing. Separate sky exposures were taken for 
each filter. The data were reduced with Standard near-IR routines 
in the IRAF UPSQIID package. The photometry is presented in 



^ HCSS is a Joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS, and SPIRE consortia 



^ Microwave Aniso tropy Dataset mapper dCantalupo et al.ll201Cll . see 
iPoditsch et"än (120 l(j ) for its implementation in HCSS) 
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Appendix|B] The K^ and H2 Images are shown in Figs.[3]and|4l 
respectively. 



3. Data analysis 

3.1. Herschel 70 /im and 160 /im 

A cursory inspection of the Herschel data (Figs.[2] and O im- 
mediately reveals a flux decrement in both PACS maps, which 
closely follows the morphology of the putative dark globule in 
NGC 1999. Following Stutz et al. (2009b), we determine the col- 
umn density and mass of the putative globule in the foreground 
of NGC 1999, assuming that the observed 70 and 160yum decre- 
ment is due to extinction. The optical depth is given by 



Table 1. Far-infrared optical depth, mass, and column density 
estimates 



T=-1.0.1n[(/ + /BG)/(/0+/BG)], 



(1) 



where / is the shadow flux level, /o is the intrinsic unabsorbed 
flux level, and /bg is the background contribution to the flux 
zero-level in the Images. We estimate the shadow flux level / 
as the mean pixel value in a ~7'' radius region centered on the 
darkest part of NGC 1999; the unobscured flux level /o is esti- 
mated as the median pixel value in a half-annulus with inner and 
outer radii of 15'' and 30'', with a southwest orientation chosen 
to avoid the bright V 380 Ori flux. The / and /o Image regions 
are over-plotted on the PACS maps in Fig.O The /bq flux level 
is not present in the Herschel maps because constant flux level 
pedestals are removed by the data reduction. We estimate the 
70yum BG flux value by interpolating the IRAS 60 and lOOyum 
Images; the 160yum BG flux values were obtained from the ISO 
Serendipity Survey in the near vicinity of NGC 1999. In both 
cases we take 50% of the total measured flux as our estimate 
for /bg as we are only interested in the missing flux originat- 
ing from hehind NGC 1999, and not the foreground component. 
While this assumption is crude, the /bg values are low compared 
to the / and /o levels and do not have a large eff'ect on our anal- 
ysis. 

Following Eq.[T] we calculate the mean optical depth per 
pixel. Table[Tlpresents the results for the data reduced with both 
methods (which we find to agree well, indicating that both re- 
cover extended emission equally well). From the resulting op- 
tical depth we calculate the column density and mass (in a 7" 
radius aperture) required to cause this flux-decrement at the 
two PACS w avelengths; these are p r esente d in Table[T] We use 
the extreme lOssenkopf & Henningl dl 9941) model dust opaci- 
ties for grains with thick ice mantles: /cyo = 541 cm^g"^ and 
/C160 = 55.7cm^g"\ a gas-t o-dust mass-r atio equal to 100, 
and a distance of 420pc (e.g.. iMenten et al.1 '2007). The result- 
ing column densities and masses for the two wavelength bands 
are clearly inconsistent. For example, we derive that masses of 
0.1 Mq and 2.5 M© are needed to account for the observed flux 
decrement in the the 70 and 160yum data, respectively. 

3.2. Ground-based follow-up 

Neither the APEX 350 nor 870//m maps show an emission fea- 
ture at the location of the Herschel flux decrement. Instead, 
the submillimeter emission bears a strong resemblance to the 
160yum emission, suggesting that the 160, 350, and 870yum maps 
are all tracing the morphology of the dust emission (Fig. [3]). This 
morphology is suggestive of a difl'use ring of material, not a hot 
region absorbed by a foreground cold cloud. Assuming optically 



Ä 


Jy/B'' 


fol fßGl T Mass 
Jy/D" Jy/D" Mo 


N(H+He) 
cm-2 


70 
160 


0.002 
0.003 


photProject (method 1) 
0.003 0.001 0.23 0.1 
0.007 0.002 0.48 2.5 


1.9 • 10^2 
3.8 • 1023 


70 
160 


0.002 
0.004 


MADmap (method 2) 
0.003 0.001 0.2 0.1 
0.008 0.002 0.45 2.4 


1.6-1022 
3.5 • 1023 



^ Derived in a ~ 7" radius region centered on R.A. = 5^36'" 24.4^ 
Decl. = -06°42'56.4", indicated in Fig. [3]as the inner most circular 
overlay in the Herschel PACS Images. 

^ Derived in a ~ 15" to ~ 30" radius half-annulus, excluding bright 
areas associated with V 380 Ori, indicated in Fig.[3]as well. 

'^ Estimated background flux contribution to the PACS maps, esti- 
mated as 50% of the measurements from IRAS data (interpolated to 
70yL/m) and the 160yL/m ISO Serendipity Survey near NGC 1999. 



thin dust emission, total gas-hdust masses can be obtained from 
submillimeter fluxes as 



Mg 
Mg+d = -rf- 



Syd^ 



Mb By(Tj))Ky^J) ' 



(2) 



where ^ is the gas-to-dust mass ratio, Sy the flux density, d the 
distance, By(Tj)) the Planck function, and /Cy,D the dust opacity. 

With an rms noise level of ~13 mJy in our LABOCA map, 
we should be able to detect (3 er) a 40 mJy point source, cor- 
responding to a mass of 0.1 3 Mp, assuming Td = 10 K and 
lOssenkopf & Henningl dl 9941) opacities for thick ice mantle 
grains. Although the 870 //m sensitivity limits are uncertain 
due to the surrounding highly structured low surface-brightness 
emission, we conclude that a mass around 0.1 M© should have 
been detected in the LABOCA map. More interestingly, the rms 
noise on the SABOCA map is around 17mJy, implying a 3cr 
point-source detection limit of ~50mJy. Due to the steeply ris- 
ing Planck function and dust opacity, this limit corresponds to a 
much lower mass detection limit of only 2.4 • 10"^ M©. For com- 
parison, the mass necessary to produce an extinction shadow as 
derived from the 70 ^m map (on the order of 0. 1 M©) would pro- 
duce a 350yum source with a flux of -200 mJy, i.e., a > lOcr de- 
tection. This result is in stark contradiction with the dense glob- 
ule/shadow Interpretation of the absorption features observed in 
the Herschel data. 

The PANIC Images show five faint stars around V 380 Ori, 
marked in Fig. [3] with their source designation as listed in Table 
IB. 11 where we also give their (H - K^) colors and estimated ex- 
tinction. Stars 3 and 4 are within the dark patch and have mod- 
erately red colors, implying a mean extinction of Ay ~ 10. In 
comparison, stars 1, 2, and 5 outside the dark patch are slightly 
redder and have a mean extinction of Ay ~ 12. None of the 
five Stars are apparent in the HST Images, indicating that they 
are background stars (Fig.[T]). The H - K^ colors show that the 
extinction toward the dark patch is much less than the 100 Ay 
needed to attenuate the 160yum emission and that the extinction 
through the dark patch is slightly less than that through the sur- 
rounding region. 

4. Discussion 

Optical Images of NGC 1999 show a dark patch suggestive of 
a small, dense globule obscuring the NGC 1999 reflection neb- 
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PANIC Ks 2,2 tirc\ 





Fig. 3. 1:75x1:75 Images of the region around NGC 1999 (including V380Ori) centered on R.A. = 5^36^24.4^ Decl. = 
-06°42'56.4'' (J2000). The white circle (14'' diameter) on the PANIC near-IR Ks Image indicates the NGC 1999 region; the Image 
is annotated with the source designations for five faint stars in the field as listed in Table IB. ll where we give their (H - K^) colors. 
The PACS 70yum and 160yum maps both show NGC 1999 observed as a flux deficit; the overlays represent the pixel masks used for 
the analysis of the NGC 1999 feature. The right-most two panels show the SABOCA and LABOCA maps. 
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Fig. 4. Left: PANIC near-IR H2 narrow-band Image of the 
NGC 1999; right: the same image with a scaled ^s continuum 
image subtracted. 



ula. Surprisingly, our 70 and 160yum Images clearly show a dark 
patch against the bright nebula with a strikingly similar mor- 
phology to that in visible light Images. These Herschel mea- 
surements, combined with subsequent ground based data, lead 
to the conclusion that the dark feature is a hole in the NGC 1999 
nebula. First, we find that the masses needed to cause the 70 
and 160yum dark patch — 0.1 and 2.5 M©, respectively — are 
inconsistent. Furthermore, the globule is not detected in emis- 
sion at 350 and 870yum: the SABOCA obervations place an Up- 
per limit of 2.4 • 10"^ M© for a temperature of 10 K. This upper 
limit is far below the amount of mass needed to cause the ob- 
scuration in the PACS data. Finally, near-IR observations with 
PANIC detect background stars toward the dark patch; the H-K^ 
colors of these stars are slightly bluer than stars detected out- 
side the globule, suggesting a lower extinction toward the dark 
patch. Furthermore, the extinctions of the background stars are 
less than that required to absorb the 160yum flux. Taken together, 
these Observation show that the dark patch is not a globule, but 
instead a cavity in the nebula. 

The presence of a well delineated cavity of the size of 
-10,000 AU deserves some attention. With the typical turbulent 
velocities on the order of a few km/s in clouds, such a cav- 
ity should be filled on timescales of at most a few 1 0,000 yrs 
and quickly disappear. The PANIC H2 narrow band data (Fig.© 
deliver important hints about the possible origin and pecu- 
liar shape of the cavity. They reveal a previously unknown, 
faint H2 bo w shock enveloping the SMZ6 0/HH 148 compact 
knots (.Stanke et al.ll2002l: Icörcoran & Ravl 11995.) . constituting 



the clearest evidence so far for a collimated flow running north- 
east to Southwest through the cavity. This flow, which likely orig- 
inales in the V 380 Ori multiple System, could possibly excavate 
the Southern part of the cavity. SMZ 6-8 in the southeastern cor- 
ner of the PANIC image resembles a small bow shock in a flow 
Coming from the north west; together with HH 35, located north- 
west of V 380 Ori (Fig.[T]), it indicates a second, northwest to 
southeast oriented flow, which could be responsible for digging 
the north western lobe of the cavity. 

We note a similar flux depression in the 160yum Herschel 
image southwest of the protostar HOPS 166, which drives the 
HH 147 outflow (ICorcoran & Ravl 1 19951) . Optical Images (e.g., 
Fig-ffl) show a circular reflection nebula marking the HH 147 
outflow cavity, coinciding with the far-IR flux depression. The 
NGC 1999 dark patch may therefore only be a somewhat pe- 
culiar example of a cavity carved in the ambient medium by 
an outflow, rendered particularly visible by the Illumination and 
heating of the cavity walls by V 380 Ori. Sensitive far-IR maps 
taken with Herschel may therefore provide a new tool to assess 
the importance of outflow feedback on cloud cores. 
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Appendix A: SABOCA and LABOCA submillimeter 
maps 

Appendix B: PANIC H and ^s photometry 

Five faint stars are visible in the PANIC H and ^s Im- 
ages. Their positions and photometry are listed in Table IB. II 
Photometrie calibration was performed with Two Micron All 
Sky Survey (2MASS|l) photometry of comparison fields con- 
taining 20 2MASS sources and verified with a 2MASS star in 
the science field. Photometry was conducted with the IRAF task 
apphot. 

The measured {H - Kg) color was used to estimate the ex- 
tinction Ay toward the stars. We assumed an intrinsic {H - Kg) 
color of 0. 17, derived as the median (H-Ks) colors of stars in an 
unreddened control field nea r Orion (Megeath et al., in prep). We 
adopted the reddening law of llndebetouw et al] (l2QQ5l) to convert 
(H-Ks) excess into ^-band extinction Ak, from which we ob- 
tained the optical extinction as Ay = 8 • A^, corresponding to a 
total to selective extinction Ry between the diff'use ISM value of 
3.1 and the higher valu es (up to 5) found for dense clouds (e.g. 
ICampeggio et al .1120071) . 



^ The Two Micron All Sky Survey is a Joint project of the 
University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 
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Fig. A.l. Füll submillimeter maps. Left: SABOCA, displayed at the original resolution of 778; right: LABOCA. Protostars in the 
field are labeled with their HOPS designations. For both maps a logarithmic color scale is used. 



Table B.l. PANIC NGC 1999 near-IR colors of background 
Stars 



Source ^ 


RA (J2000) 


Decl. (J2000) 


Ks 


cr{Ks) 


(H-Ks) 


o-(H-Ks) 




(hms) 


r ' ") 


mag 


mag 


mag 


mag 


1 


05 36 27.54 


-06 43 22.29 


16.08 


O.Ol 


0.83 


O.Ol 


2 


05 36 26.19 


-06 42 46.23 


15.68 


O.Ol 


0.95 


0.02 


3 


05 36 24.78 


-06 43 06.72 


18.36 


0.10 


0.67 


0.14 


4 


05 36 24.09 


-06 42 51.12 


16.62 


O.Ol 


0.84 


0.02 


5 


05 36 23.48 


-06 42 51.86 


16.71 


0.02 


0.95 


0.02 



" See Fig.[3]for source designation. 



